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Solar  Radiation  Alert  System 


INTRODUCTION 

Aircraft  passengers  and  crew  are  constantly  exposed 
to  ionizing  radiation  at  higher  dose  rates  than  normally 
received  by  the  general  population.  The  principal  ion¬ 
izing  radiation  is  galactic  cosmic  radiation,  a  main  source 
of  which  is  believed  to  be  supernovae  (exploding  stars). 
Occasionally,  a  disturbance  in  the  sun  (solar  flare,  coronal 
mass  ejection)  leads  to  a  large  flux  of  solar  protons  with 
sufficient  energy  to  penetrate  the  earth’s  magnetic  field, 
enter  the  atmosphere,  and  increase  ionizing  radiation 
levels  at  aircraft  flight  altitudes.  For  a  primer  on  aircrew 
exposure  to  ionizing  radiation,  see  report  by  Friedberg 
and  Copeland  (1). 

A  system  (Figure  1)  has  been  developed  for  the  con¬ 
tinuous  evaluation  of  proton  flux  measurements  made 
by  instruments  on  Geosynchronous  Operational  Envi¬ 
ronmental  Satellites  (GOES).  At  present,  the  primary 
data  source  is  GOES- 1 1  while  GOES- 10  measurements 
are  used  if  they  are  more  current  or  if  GOES- 1 1  is  not 
operating  properly.  If  the  measurements  indicate  the 
likelihood  of  a  substantial  elevation  of  ionizing  radia¬ 
tion  levels  at  aircraft  flight  altitudes  at  high  latitudes, 
the  Civil  Aerospace  Medical  Institute  (CAMI)  issues  a 
Solar  Radiation  Alert  (SRA)  to  the  National  Oceanic 
and  Atmospheric  Administration  Weather  Wire  Service 
(NWWS)  (2).  Information  given  with  the  alert  includes 
estimated  dose  rates  in  microsieverts  per  hour  (pSv/h), 
based  on  5-minute  averages,  at  30,000,  40,000,  50,000, 
60,000,  and  70,000  ft  (relative  to  mean  sea  level).  The 
dose  rates  are  updated  at  5-minute  intervals  for  the  dura¬ 
tion  of  the  SRA. 

The  specific  criteria  for  issuance  of  an  SRA  are  that 
the  estimated  effective  dose  rates  induced  by  solar 
protons  at  any  altitude  from  30,000  to  70,000  ft  in 
the  earth's  atmosphere  are  equal  to  or  greater  than  20 
microsieverts  per  hour  for  each  of  three  consecutive  5- 
minute  periods.  The  three  5-minute  periods  need  not 
be  at  the  same  altitude.  Because  of  the  potential  for  a 
delay  of  3  to  20  minutes  in  communication  processes 
and  the  three  5-minute  periods,  notification  of  the  start 
of  SRA  conditions  may  be  18  to  35  minutes  after  the 
start  of  the  first  5 -minute  period.  Possible  effects  on  the 
dose  rate  caused  by  ions  other  than  protons  (e.g.,  alpha 
particles)  in  the  primary  solar  particle  flux,  shielding 
by  the  earth's  magnetic  field,  and  shielding  by  aircraft 
structure  are  not  taken  into  account  in  the  calcula¬ 


tions.  The  SRA  is  canceled  when  the  average  effective 
dose  rate  at  all  altitudes  from  30,000  to  70,000  ft  is 
less  than  20  microsieverts  per  hour  for  six  consecutive 
5-minute  periods. 

METHOD  USED  TO  ESTIMATE 

eefecttve  dose  rates  at  flight 
altitudes 

Basis  of  SRA  System 

The  Primary  Solar  Proton  Spectrum.  Observa¬ 
tions  of  fluxes  of  high-energy  protons  used  to  estimate 
solar  radiation  dose  rates  in  the  earth’s  atmosphere  are 
obtained  from  GOES  proton  flux  measurements.  These 
measurements  are  recorded  in  several  “channels”:  P4  to 
P7  (15  to  900  MeV  energy,  168  to  1581  MV  rigidity), 
PIO  (510  to  700  MeV,  1103  to  1343  MV),  and  Pll 
(greater  than  700  MeV,  greater  than  1343  MV).  Protons 
in  channels  P4  through  PIO  are  in  terms  of  differential 
flux  (particles/cm^/steradian/second/MeV).  Protons  in 
channel  Pll  are  in  terms  of  integral  flux  (particles/cm^/ 
steradian/second). 

The  proton  fluxes  currently  available  from  the  Space 
Environment  Genter  (SEG)  are  derived  from  satellite- 
instrument  count  rates  using  conversion  factors  ap¬ 
propriate  for  GOES-7  instruments  (3).  In  Steps  1  to  3, 
below,  these  fluxes  are  converted  to  count  rates  and  then 
reconverted  to  flux  using  the  conversion  factors  specifi¬ 
cally  intended  for  GOES-10  and  GOES-11.  Ghannels 
P8  and  P9  are  excluded  from  the  analysis  because  their 
correction  algorithms  are  more  complex  than  those  of  the 
other  channels,  and  excluding  them  does  not  significantly 
affect  the  calculated  dose  rates. 

The  data  are  available  in  near  real-time  from  the  file 
transfer  protocol  (ftp)  server  at  the  SEG  and  are  also 
on  the  Internet  (4).  Protons  flux  observations  are  from 
solar  plus  galactic  cosmic  radiation.  The  contribution  of 
galactic  protons  is  subtracted  from  the  total  as  part  of  the 
correction  process  described  below.  The  initial  5-minute 
average  galactic  background  count  rate  {B,  counts/ second) 
for  a  channel  is  the  average  of  the  preceding  12  hours 
of  data  (144  5-minute  periods)  collected  when  no  solar 
particle  event  or  Forbush  decrease  (sudden  reduction 
in  galactic  background  radiation  as  a  result  of  scatter¬ 
ing  of  low-energy  galactic  particles  by  magnetic  fields 
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carried  by  solar  particles)  is  occurring.  Once  the  system 
is  operating,  the  background  is  changed  as  described  in 
Step  10,  below. 

In  each  channel,  solar  protons  of  sufficiently  high  en¬ 
ergy  to  cause  an  increase  in  ionizing  radiation  levels  in  the 
atmosphere  have  a  flux  spectrum  (f)  across  the  channel’s 
energy  acceptance  band  that  can  reasonably  be  expressed 
as  a  power  law  function  of  rigidity  (R): 

f(R)  =  aR^  (Eq.  1) 

where  the  proton  spectrum  is  characterized  in  terms  of 
time-dependent  constants,  intensity  (a)  and  spectral 
hardness  index  (y)  (5,6).  Six  spectra  in  the  form  of  Eq.  1 , 
each  spanning  a  different  rigidity  range,  are  created  in 
Steps  1  to  7  (below)  using  GOES  data. 

These  six  spectra  are  combined  in  Step  8  to  form  a  com¬ 
plete  spectrum  used  to  calculate  effective  doses  in  Step  9.  A 
piecewise-continuous  spectrum  is  needed  because  during 
solar  proton  events,  y  increases  at  higher  rigidities. 

Radiation  Transport.  Monte  Carlo  radiation  transport 
code  MCNPX2.4.0  (7)  was  used  to  calculate  the  fluences 
(particles/cm^)  of  secondary  particles  (protons,  neutrons, 
pions,  kaons,  photons,  electrons,  and  muons)  in  selected 
energy  ranges  for  speciflc  altitudes,  which  were  generated 
by  mono-energetic  isotropic  fluences  of  protons  at  an 
altitude  of  100  km.  Solar  proton  fluences  at  100  km  were 
assumed  to  be  isotropic  because  there  is  no  directional 
information  in  the  GOES  data.  The  atmospheric  pres¬ 
sure  density  at  100  km  is  0.0003  g/cm^.  Using  a  higher 
altitude  to  more  closely  approximate  the  free-space  envi¬ 
ronment  of  the  satellite  would  not  significantly  improve 
dose  estimates.  The  earth  was  modeled  as  a  sphere  of 
liquid  water  with  a  radius  of  6371  km  and  a  density  of 
Ig/cm^.  The  atmosphere  was  modeled  as  a  series  of  100 
spherically  symmetric  1-km  thick  layers  surrounding  the 
earth.  The  composition  (a  gaseous  mixture  containing 
nitrogen,  oxygen,  argon,  and  carbon  atoms)  and  density 
of  each  layer  were  assumed  to  be  that  reported  for  the 
middle  of  the  layer  (8),  resulting  in  an  atmosphere  with 
a  vertical  depth  of  1035.08  g/cm^.  Altitudes  in  Table  1 
were  calculated  from  the  atmospheric  depth.  Geomagnetic 
effects  were  neglected. 

Fluence-to-Eflfective  Dose  Conversion.  With  one 
exception,  fluence-to-effective-dose  conversion  coef¬ 
ficients  reported  by  Pellicioni  (9)  were  used  to  convert 
the  secondary  particle  fluences  generated  by  the  transport 
code  to  effective  doses.  The  reported  coefficients  incor¬ 
porate  radiation-weighting  factors  and  tissue -weighting 
factors  recommended  by  the  International  Gommission 
on  Radiological  Protection.  The  one  exception  is  the 


radiation-weighting  factor  for  protons,  which  was  changed 
from  five  to  two,  as  recommended  by  the  National 
Gouncil  on  Radiation  Protection  and  Measurements 
(10).  The  effective  dose  from  kaons  was  not  calculated 
because  its  contribution  to  the  total  dose  would  have 
been  negligible.  Table  1  shows,  at  selected  altitudes,  the 
dose  rate  per  unit  flux  of  primary  solar  protons  of  selected 
rigidities.  The  transport  code  automatically  converts  a 
fluence  below  10^^^  particles/cm^  to  zero  fluence,  which 
would  lead  to  a  zero  dose  rate  per  unit  flux  (Table  1).  To 
avoid  this,  any  zero  dose  rate  per  unit  flux  was  changed 
from  zero  to  1x10 Table  1,  as  presented,  is  applicable 
to  any  solar  particle  event. 

Calculating  Effective  Dose  Rate  at  Specific  Altitudes 

Step  1.  Gonvert  the  incorrect  fluxes  obtained  from 
the  SEG  to  instrument  count  rates  (Q,  by  multiplying 
the  flux  (^)  for  each  channel  used  (pp^,  (p^^,  (p^^, 

Ppjg,  and  Ppjj)  by  the  appropriate  conversion  factor  {k) 
from  Table  2; 

C  =  (pk  (Eq.  2) 

GOES- 10  and  GOES- 1 1  have  identical  instruments  and 
therefore  use  the  same  conversion  factors. 

Step  2.  Subtract  the  galactic  background  count  rate 
(B)  from  the  total  count  rate  (C,  galactic  plus  solar)  to 
obtain  the  solar-induced  count  rate  {C'): 

C=  C-B  (Eq.  3) 

If  C  is  negative,  assign  it  a  value  of  1%  of  B.  Galcula- 
tion  of  B  is  described  earlier,  in  the  section  entitled  “The 
Primary  Solar  Proton  Spectrum.” 

Step  3.  Divide  the  background  corrected  solar  count 
rate  {C)  by  the  appropriate  conversion  factor  (k')  from 
Table  2  to  obtain  the  differential  flux  (f(R)): 

f(R)  =  C/k'  (Eq.  4) 

Step  4.  Galculate  a  preliminary  spectral  hardness  index 
(y)  for  each  pair  (i,j)  of  channels  (P4  and  P5,  P5  and  P6, 
P6and  P7,P7  andP10,P10andPll)  using  Eq.  5, which 
was  derived  from  Eq.  1 ; 

Tpi,pj  1  /p  \  1  /p  y  ^ 

Variables  f(R)p^  through  f(R)p^^  are  differential  proton 
fluxes  derived  from  values  in  channels  P4  through  P 1 1 
(Steps  1  to  3).  The  characteristic  rigidities  of  the  protons 
in  these  channels  are  i?p^  through  i^p^^  from  Table  2. 
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Step  5.  Further  correct  the  differential  fluxes  (f(R)j,.J, 
using  spectral  hardness  index  correction  factors  (s)  in 
Table  3: 

(Eqs) 

Use  the  p.  values  from  Step  4  to  And  the  appropriate 
s  values.  Thus,  use  y^,^  p^  to  And  the  s  value  to  correct 
use  yp^pg  to  And  the  s  value  to  comet  f(R)^^,  etc. 
This  selection  procedure  weights  the  spectral  response 
towards  the  higher  energies. does  not  require  this 
correction. 

Step  6.  Recalculate  y  values  with  Eq.  5,  using  fluxes 
fiR)^^  and/f^jp.  that  were  corrected  in  Step  5.  If  y  is 
negative,  assign  it  the  value  0.  If  y  equals  I,  assign  it  the 
value  I.OOOOOI. 

Step  7.  Calculate  a  for  each  channel  with  Eq.  7  using 
the  characteristic  rigidities  from  Table  2,  the  differential 
fluxes/(7?j^^_^_^__^^^j  from  Step  5 ,  and  the  recalculated  y  values 
from  Step  6: 

“  7) 

with/f7?Jp^  use  pp  with  use  etc.  Fot  chan- 
nel  PI  1,  use  the  a  found  for  PIO.  Convert  channel  PI  1 
count  rate  ( C)  to  integral  flux  (Tp^j),  irrespective  of  y,  by 
dividing  C  by  its  geometric  factor  0.73  (3).  Use  Eq.  8  to 
And  the  y  that  best  describes  .Fp^j.  Start  with  y  equal  to 
ypjppjj  and,  if  necessary,  decrease  y  in  steps  of  0.000001 
until  the  relationship  is  just  satisfied, 
c  32545MF 

Fpi,  <  a  R  ^ dR  (Eq.  8) 

The  upper  and  lower  limits  in  Eq.  8  are  the  rigidity  bounds 
set  for  channel  P 1 1 .  The  upper  limit  is  the  upper  rigidity 
limit  of  the  Monte  Carlo  calculations,  32545  MV  (Table 
1),  and  the  lower  limit  is  1343  MV,  the  lowest  rigidity 
of  data  in  P 1 1 . 

Step  8.  Construct  a  piecewise-continuo us  description 
of  the  solar  proton  spectrum  (nf(R))  for  primary  solar 
protons  with  rigidities  above  137  MV  (10  MeV)  based 
on  results  of  Steps  1  to  7,  with  each  piece  having  the 
form  of  Eq.  1 .  With  the  exception  of  the  upper  rigidity 
limit  of  channel  P 1 1 ,  the  upper  and  lower  rigidity  limits 
used  with  the  parameters  calculated  for  each  channel  are 
shown  in  Table  4.  The  upper  rigidity  limit  for  channel 
P 1 1  (and  thus  for  \ii(R))  is  R^^  and  is  found  by  replacing 
the  upper  rigidity  limit  in  Eq.  8  with  1344  MV,  then 
increasing  the  upper  limit  in  steps  of  1  MV  until  the 
relationship  is  just  satisfied. 


Step  9.  For  a  specific  altitude  (H),  the  effective  dose 
rate  per  unit  flux  of  primary  solar  protons  {^s.(R,H))  is 
generated  from  the  data  in  Table  1 .  For  intermediate  R 
values,  A(R,H)  is  generated  using  linear  interpolation  of 
Log(R)  and  LogiA(R,H)).  Primary  solar  protons  in  bins 
with  rigidities  less  than  183.5  MV  (17.78  MeV),  the 
lowest  rigidity  in  Table  1,  are  assigned  the  same  A(R,H) 
as  primary  solar  protons  of  183.5  MV.  The  effective 
dose  rate  (E(H),  microsieverts  per  hour)  at  a  particular 
altitude  is  the  integral 

E{H)=  \  y/{R)A{R,H)dR  (Eq.  9) 

where  R  .  is  the  lower  limit  of  channel  P4  (Table  4) 
and  R  is  from  Step  8.  The  numerical  integration  is 
performed  using  1  MV-wide  flux  bins.  For  each  bin,  the 
A(R,H)  value  used  is  the  mean  of  the  A(R,H)  values  at 
the  extreme  rigidities  of  the  bin. 

Dose  rates  at  altitudes  between  those  listed  in  Table  1 
are  estimated  using  a  cubic  spline  interpolation. 

Step  10.  If  C  from  Step  2  does  not  exceed  twice  B  in 
any  channel,  C  is  used  to  update  B  in  each  channel: 

B=  (143B+C)/144  (Eq.  10) 

Thus,  the  background  count  rates  will  always  be  the 
averages  of  the  last  12  uneventful  hours. 

To  examine  solar  proton  events  using  proton  flux 
data  from  satellites  other  than  the  current  GOES  series 
(GOES-8,  GOES-9,  GOES-10,  GOES-1 1,  GOES-12), 

minor  modifications  are  necessary;  any  satellite  data  set 
may  be  used  that  allows  the  reasonable  construction  of 
a  solar  proton  spectrum  spanning  the  rigidity  range  137 
MV  to  above  1343  MV. 

RESULTS  AND  DISCUSSION 

In  estimating  dose  rates  for  the  SRA  system,  only 
solar  protons  are  used  because  they  are  the  dominant 
contributor  to  the  ionizing  radiation  from  the  sun  at 
aircraft  flight  altitudes  and  data  are  available  in  near 
real-time.  Other  radiation  information,  such  as  GOES 
alpha  particle  measurements  and  data  from  the  worldwide 
network  of  neutron  monitors,  would  enable  more  com¬ 
plete  calculations,  but  not  with  the  immediacy  required 
of  an  alert  system. 

During  the  early  stages  of  a  solar  proton  event  (solar 
particle  event  where  the  proton  flux  with  energies  equal 
to  or  greater  than  10  MeV  exceeds  10  particles/second/ 
steradian/cm^  (1 1)),  solar  radiation  in  the  vicinity  of  the 
earth  is  usually  somewhat  anisotropic,  as  indicated  by 
ground-level  neutron  monitor  data  (12).  Flowever,  an¬ 
isotropy  is  not  particularly  evident  in  GOES  data.  The 
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15  largest  solar  proton  events  from  1996  through  2003 
for  which  data  were  available  from  two  or  more  satellites 
were  examined  to  compare  event  onset  times  and  flux 
measurements.  Apparent  event  onset  times  at  the  satellites 
were  the  same  (within  5  minutes).  The  largest  differences 
in  fluxes  measured  by  the  satellites,  typically  20%,  were 
observed  at  energies  greater  than  350  MeV  during  the  first 
few  5-minute  periods  of  an  event.  As  the  event  proceeded, 
differences  decreased.  Longitudinal  separation  of  the  most 
widely  separated  satellites  was  about  60  degrees.  The  low 
anisotropy  effects  observed  at  the  satellites  indicate  that 
doses  generated  by  the  system  are  reasonable  first  estimates 
at  aircraft  flight  altitudes  at  polar  latitudes  (i.e.,  where 
the  cutoff  rigidity  is  close  to  0  MV). 

Figure  2  shows  effective  dose  rates  from  solar  protons 
in  the  earth’s  atmosphere  during  the  solar  proton  event 
of  20  January  2005.  The  calculations  do  not  take  into 
account  geomagnetic  shielding  (i.e.,  shielding  by  the 
earth’s  magnetic  held)  and  shielding  by  aircraft  structure. 
Inclusion  of  aluminum  shielding  afforded  aircraft  occu¬ 
pants  in  a  B-747  (using  a  skin  thickness  of  0.6  g/ cm^) ,  for 
example,  would  not  change  the  effective  dose  estimates  at 
40,000  ft  by  more  than  a  few  percent  (depending  on  the 
specific  energy  spectrum,  the  effective  dose  could  increase 
or  decrease).  The  event  started  at  06:50  Universal  Time 
(UT),  and  by  06:55  UT  dose  rates  at  60,000,  40,000, 
and  30,000  ft  had  reached  maximum  values  of:  140,  55, 
and  2 1  microsieverts  per  hour,  respectively.  The  dose  rates 
then  dropped  rapidly:  at  07:55  UT,  they  were  17,  4,  and 
1  microsieverts  per  hour  at  the  three  altitudes. 

The  SRA  system  provides  timely,  quantitative  in¬ 
formation  for  the  aviation  community  on  potentially 
hazardous  solar  radiation  levels  at  aircraft  flight  altitudes 
in  high-latitude  regions.  This  information  enables  users 
to  decide  whether  or  not  reducing  aircraft  flight  altitudes 
is  appropriate. 

Future  work  includes:  (a)  improve  estimates  of  dose 
rates  in  Table  1  by  improving  Monte  Carlo  calculations, 
(b)  improve  estimates  of  radiation  levels  at  neutron  moni¬ 
tor  locations  by  including  geomagnetic  and  anisotropy 
effects  on  dose  rates,  (c)  possible  incorporation  of  dose 
rates  from  alpha  particles  and  any  other  significant  heavier 
ions  in  the  solar  particle  flux,  and  (d)  calculation  of  dose 
rates,  flight  doses,  and  neutron  monitor  counting  rates 
during  past  solar  proton  events. 
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FIGURES  AND  TABLES 


Figure  1.  Solar  radiation  alert  system:  (1)  Occasionally,  a  disturbance  in  the  sun  (solar  flare,  coronal  mass 
ejection)  leads  to  a  large  flux  of  high-energy  particles  in  the  vicinity  of  the  earth.  (2)  Instruments  on  a  GOES 
satellite  continuously  measure  the  radiation  and  the  information  is  transmitted  to  NOAA.  (3)  From  there  it 
is  sent  to  CAMI.  A  computer  at  CAMI  analyzes  the  measurements.  (4)  If  the  measurements  indicate  the 
likelihood  of  a  substantial  elevation  of  ionizing  radiation  levels  at  aircraft  flight  altitudes,  a  Solar  Radiation 
Alert  is  issued  to  the  NOAA  Weather  Wire  Service  within  10  minutes.  (5)  NOAA  Weather  Wire  Service 
subscribers  are  provided  estimated  effective  dose  rates  for  30,000,  40,000,  50,000,  60,000,  and  70,000  ft. 
This  information  is  updated  at  5-minute  intervals  for  the  duration  of  the  alert. 
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Figure  2.  Effective  dose  rates  from  solar  proton-induced  ionizing  radiation  at 
selected  altitudes  on  20  January  2005  (based  on  GOES-11  data,  geomagnetic 
effects  neglected). 
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Table  1.  Effective  dose  rate  per  unit  flux  of  primary  solar  protons,  as  related  to  incident  rigidity 
and  energy  at  selected  altitudes  in  the  atmosphere  relative  to  mean  sea  level. 


Dose  rate  per  unit  flux^  A(R,H),_  microsieverts  per 
hour/(protons/cmVsteradian/second)  at  altitude  (H) 

Rigidity  (R),  MV 

Energy,  MeV 

H=  -50  ft 

77=9722  ft 

77=19,588  ft 

77=29,378  ft 

183.5 

17.78 

1.000x10'** 

1.000x10'** 

1.000x10'** 

2.901x10''* 

245.6 

31.62 

1.000x10'** 

1.000x10'** 

1.404x10'*** 

1.256x10'* 

329.7 

56.23 

1.000x10'** 

1.000x10'** 

2.978x10'^ 

2.063x10'" 

444.6 

100 

1.000x10'** 

9.739x10'^ 

2.251x10'" 

7.796x10'® 

604.4 

177.8 

9.040x10'* 

1.919x10'" 

7.713x10'® 

1.271x10''* 

832.7 

316.2 

5.281x10'" 

2.489x10'* 

2.935x10'^ 

1.960x10'* 

1171. 

562.3 

7.556x10'® 

2.383x10'^ 

2.267x10'* 

1.207x10'* 

1696. 

1000 

4.294x10'* 

1.122x10'* 

9.903x10'* 

5.130x10'* 

2549. 

1778 

2.173x10''* 

4.428x10'* 

3.625x10'* 

0.1469 

3991. 

3162 

7.719x10'^ 

1.495x10'* 

9.106x10'* 

0.3009 

4848. 

4000 

1.122x10'* 

1.915x10'* 

0.1147 

0.3651 

6495. 

5623 

1.885x10'* 

2.356x10'* 

0.1319 

0.4335 

10898. 

10000 

7.433x10'* 

4.966x10'* 

0.2413 

0.7315 

18695. 

17780 

2.202x10'* 

0.1064 

0.4464 

1.277 

20917. 

20000 

2.751x10'* 

0.1227 

0.5025 

1.410 

32545. 

31620 

6.046x10'* 

0.2178 

0.7942 

2.110 

//=3  9, 197  ft 

H=4S,9n  ft 

77=58,749  ft 

77=68,511  ft 

77=78,260  ft 

183.5 

8.348x10'* 

1.282x10'® 

7.556x10'® 

1.937x10'* 

3.215x10'* 

245.6 

9.257x10'" 

1.521x10'* 

7.390x10'* 

1.699x10'^ 

2.622x10''* 

329.7 

1.199x10'^ 

1.223x10''* 

4.873x10'^ 

1.004x10'* 

1.456x10'* 

444.6 

1.511x10''* 

7.832x10'^ 

2.407x10'* 

4.436x10'* 

6.602x10'* 

604.4 

1.046x10'* 

4.495x10'* 

1.129x10'* 

1.165x10'* 

2.446x10'* 

832.7 

8.491x10'* 

2.470x10'* 

5.079x10'* 

0.1116 

0.2307 

1171. 

4.331x10'* 

0.1298 

0.2748 

0.3850 

0.4490 

1696. 

0.1731 

0.3447 

0.4922 

0.5712 

0.5900 

2549. 

0.3629 

0.5957 

0.7517 

0.7941 

0.7719 

3991. 

0.6252 

0.9125 

1.064 

1.077 

1.004 

4848. 

0.7475 

1.079 

1.243 

1.242 

1.150 

6495. 

0.9214 

1.379 

1.621 

1.653 

1.518 

10898. 

1.484 

2.100 

2.357 

2.296 

2.075 

18695. 

2.418 

3.234 

3.451 

3.241 

2.835 

20917. 

2.638 

3.499 

3.711 

3.450 

3.006 

32545. 

3.792 

4.795 

4.884 

4.425 

3.762 

®  Geopotential  altitude  in  feet  is  converted  from  atmospheric  depth  in  g/cm^  using  ref.  8. 
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Table  2.  Conversion  factors  (k,  k')  and  characteristic  rigidities  for 
measurements  in  GOES-10  and  GOES-11  proton  channels  (3) 


Channel 

Conversion 
Factor  k  ^ 

Conversion 

Factor  k' 

Characteristic 
Rigidity,  MV  ^ 

P4 

4.64 

22.25 

225.1 

P5 

15.5 

43.04 

338.2 

P6 

90. 

252.8 

563.9 

P7 

300. 

1210. 

950.0 

PIO 

162. 

175.6 

1225. 

Pll 

1565. 

1103. 

1700. 

““  k.  counts/(particles/cm^/steradian/MeV) 

'"k'.  counts/(particles/cm^/steradian/MV) 

“^When  calculating  a  spectral  hardness  index,  a  single  rigidity  value,  called  the 
characteristic  rigidity,  is  assigned  to  the  entire  channel. 


Table  3.  Spectral  hardness  index  correction  factors  (s)  for 
differential  fluxes  derived  from  measurements  in  GOES-1 0 
and  GOES-11  proton  channels  (3)  ^ 


Y 

P4 

P5 

P6 

P7 

PIO 

0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1 

0.9565 

0.9775 

0.9733 

0.9300 

0.9968 

2 

0.8777 

0.9350 

0.9233 

0.8050 

0.9904 

3 

0.7751 

0.8764 

0.8553 

0.6502 

0.9808 

4 

0.6609 

0.8060 

0.7753 

0.4913 

0.9683 

5 

0.5458 

0.7284 

0.6888 

0.3485 

0.9529 

6 

0.4380 

0.6476 

0.6011 

0.2329 

0.9348 

7 

0.3428 

0.5674 

0.5161 

0.1473 

0.9142 

8 

0.2623 

0.4904 

0.4366 

0.08856 

0.8914 

9 

0.1969 

0.4186 

0.3645 

0.05089 

0.8666 

10 

0.1453 

0.3533 

0.3006 

0.02811 

0.8400 

Intermediate  correction  factors  are  obtained  by  linear  interpolation. 


Table  4.  Lower  and  upper  rigidity  limits  associated  with  parameters  «and 
y  calculated  for  each  channel. 


P4 

P5 

P6 

P7 

PIO 

Pll 

Lower  Rigidity 

137 

338 

564 

950 

1225 

1343 

Upper  Rigidity 

338 

564 

950 

1225 

1343 

See  Step  8 
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